#ifndef COMPLEX_H
#define COMPLEX_H
/* Define complex multiplication and its conjugate */
#define  rmul(x,y)      (x.re * y.re - x.im * y.im)
#define  imul(x,y)      (x.im * y.re + x.re * y.im)
#define rcmul(x,y)      (x.re * y.re + x.im * y.im)
#define icmul(x,y)      (x.im * y.re - x.re * y.im)
#define PI  3.14159265
class Complex
{
public:
    Complex();
    float re;
    float im;
};

#endif // COMPLEX_H
